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SIMULATING  SMOKE  MOVEMENT  THROUGH  LONG  VERTICAL  SHAFTS  IN  ZONE-TYPE 

COMPARTMENT  FIRE  MODELS 


ABSTRACT 

A limitation  of  traditional  zone-type  compartment  fire  modeling  concepts  is  identified;  namely,  the 
inadequacy  of  two-layer  quasi-steady-buoyant-plume  analyses  to  simulate  the  fire-generated 
environment  in  room  configurations  with  large  height-to-span  ratios,  e.g.,  elevator  shafts  and  long, 
vertical,  ventilation  shafts  and  ducts.  A possible  means  of  removing  this  limitation  is  developed.  This 
involves  a method  of  analysis  and  associated  model  equations  that  can  be  implemented  and  used  to 
advance  zone-type  models.  The  model  equations  simulate  time-dependent  flows  in  a long,  ventilated, 
vertical  shaft/duct  with  an  arbitrary  vertical  density  distribution,  including  one  or  more  intervals  along 
the  shaft/duct  length  where  the  vertical  distribution  of  the  average  cross-section  density  may  be 
unstably  stratified,  i.e.,  density  increasing  with  increasing  elevation.  The  model  equations  are  partially 
verified  by  favorable  comparisons  between  solutions  and  previously  published  data  from  unsteady 
experiments  in  long  vertical  tubes  involving  initially  unstable  configurations:  salt-water  over  fresh- 
water and  heavy-gas  over  light-gas.  Additional  verification  of  the  proposed  equation  set  with  cold-air 
over  hot-air  systems  and  with  fire-driven  smoke  flows,  both  of  which  involve  gas-to-surface  heat  trans- 
fer, is  required  before  this  model  can  be  used  with  confidence  in  professional  practice. 


INTRODUCTION 

Zone-Type  Compartment  Fire  Modeling  Concepts 

Zone-type  computer  fire  models  are  important  practical  tools  for  simulating  mathematically  the 
dynamics  of  fire  and  smoke  spread  through  multi-room  facilities.  In  such  models,  rooms  are  assumed 
to  be  connected  to  one  another  or  to  the  outside  environment  by  vents  (i.e.,  holes)  in  bounding 
partitions.  Room-to-room  cross-vent  exchanges  of  fresh  air  and/or  smoke  are  determined  from 
momentum  considerations  (e.g.,  the  Bernoulli  equation);  these  exchanges  being  driven  by  room-to- 
room,  cross-vent,  hydrostatic  pressure-differences.  The  vertical  hydrostatic  pressure  distribution  in 
any  particular  room  of  a modeled  facility  is  calculated  from  the  time-dependent  pressure  at  the  floor 
(or  at  some  other  fixed  elevation  in  the  room)  and  from  the  floor-to-ceiling  distribution  of  density. 
Time-dependent  pressure  at  the  reference  elevation  and  density  distribution  in  each  room  of  the 
facility  ^re  included  as  variables  computed  by  the  fire  model.  Variations  in  pressures  and  flows 
throughout  the  facility  are  determined  from  the  buoyancy  of  fire-generated  layers  of  smoke  that  grow 
in  thickness  and  temperature  at  the  upper  portions  of  one  or  more  rooms  of  the  facility,  by  non- 
uniform,  wind-generated,  stagnation  pressure  distributions  that  surround  the  ("leaky")  external 
surfaces  of  the  facility,  and/or  by  stack  effects  (elevation-dependent  pressure  differences  between  the 
inside  and  outside  of  a tall  facility  generated  by  differences  of  outside  density/temperature  and 
climate-controlled  inside  density/temperatures.  Some  fire  models  also  have  a capability  for  simulating 
the  effect  on  pressures  and  flows  of  operating  ventilation  (HVAC)  system  components.  In  the  latter 
case,  the  current  state  of  zone-type  fire  modeling  technology  is  such  that  duct  pressure  distributions 
and  duct  flow  and  mixing  are  generally  calculated  or  specified  without  a full  accounting  of  the 
buoyancy  of  the  flowing  gas. 

The  fire-generated  environment  in  any  particular  room  of  a facility  is  generally  approximated  by  two 
stably  stratified  layers;  an  upper,  elevated-temperature,  and  relatively-buoyant  smoke  layer  and  a 
cooler  and  less-buoyant  lower  layer.  Each  layer  is  assumed  to  have  spatially-uniform  time-dependc.at 


properties.  Prior  to  the  time  when  smoke  is  introduced  to  a room,  either  from  a fire  contained  in 
the  room  or  by  smoke  inflow  through  a vent  or  ventilation  duct  inlet,  the  room  would  be  filled  with 
an  ambient-temperature  lower  layer,  i.e.,  the  interface  separating  the  upper  layer  from  the  lower  layer 
would  be  at  its  initial,  ceiling  elevation. 

When  a source  of  buoyant  smoke  is  introduced  to  a room  below  the  smoke  layer  it  rises  as  a buoyant 
plume  toward  the  smoke  layer  interface,  entraining  laterally  lower-layer  gases  during  its  ascent.  In 
the  current  state  of  technology  of  zone-type  models,  the  flow  dynamics  above  a source  of  buoyant 
smoke  is  modeled  mathematically  as  a purely-buoyant  quasi-steady  plume  in  a large  quiescent 
environment.  Here,  "large"  is  intended  to  indicate  that  the  volume  of  the  lower  layer  is  large 
compared  to  the  volume  containing  the  significant  plume-flow  activity.  Also,  implicit  in  the  uniform 
layer  approximation  is  the  assumption  that,  when  a plume  enters  a layer,  the  characteristic  time  of 
mixing  of  the  penetrating  plume  gases  with  the  bulk  of  the  receiving  layer  is  small  compared  to 
characteristic  times  of  the  overall  phenomena  being  simulated. 

A Limitation  of  the  TVvo-Layer  Quasi-Steady-Buoyant-Plume  Modeling  Approach 

The  modeling  strategy  which  uses  the  concepts  of  one-to-two  uniform  layers  per  room,  room-to-room 
mass  exchanges  by  vent  flows,  and  layer-to-layer  mass  exchange  by  quasi-steady  buoyant  plumes  has 
proven  to  be  very  robust.  However,  there  are  important  practical  room  configurations  and  associated 
fire  scenarios  where  these  basic  concepts  are  inadequate.  If  the  basic  concepts  are  not  applicable  to 
a particular  room  configuration,  then,  to  the  extent  that  room  plays  a significant  role  in  the  spread 
of  fire  and  smoke  throughout  the  facility  (e.g.,  the  room  is  the  connecting  flow  path  between  the 
room  of  fire  origin  and  a threatened  space),  the  inadequacy  of  the  simulation  in  that  room  can  render 
inadequate  the  entire  simulation. 

Let  L and  d,  be  the  characteristic  height  and  horizontal  span,  respectively,  of  a room  configuration. 
Then,  in  view  of  the  above  discussion,  for  rooms  with  Ljd  > > 1 (shafts  or  ducts)  it  is  not  reasonable 
to  expect  a traditional,  two-layer,  zone-type,  fire  modeling  approach  to  lead  to  a successful  simulation 
of  fire  environments  if  such  spaces  are  described  and  treated  in  the  manner  of  a standard,  two-layer, 
room  element,  with  plumes,  etc.  For  such  room  configurations,  some  of  the  previously  mentioned 
basic  modeling  assumptions  become  invalid.  Thus,  1)  as  the  plume  rises  and  spreads,  its  volume 
eventually  becomes  significant  and  it  starts  to  fill  a large  fraction  of  the  section  of  the  shaft/duct,  2) 
it  is  not  reasonable  to  expect  that  characteristic  times  of  mixing  in  the  shaft  will  generally  be  small 
compared  to  characteristic  times  of  interest,  and  3)  there  is  no  basis  for  support  for  a uniform,  two- 
layer  approximation  to  the  density/temperature  distribution.  An  illustration  of  a generic  problematic 
facility  and  room  configuration  is  presented  in  Figure  1.  In  the  facility  shown,  traditional  modeling 
concepts  would  typically  be  applicable  in  the  small- lyj  room  of  fire  origin,  the  upper  floors,  etc.,  but 
they  would  not  be  applicable  in  the  large- I/J  shaft. 

A Strategy  for  Modeling  Flow  Through  Shaft-Like  Spaces  in  2^ne-'iype  Compartment  Fire  Models; 
Objective  of  This  Work 

For  the  traditional  two-layer,  etc.,  zone-type  modeling  approach  to  be  valid,  large- I/J  spaces  in  a 
multi-room  facility  need  to  be  treated  as  a special  class  of  room  configuration.  Fire-generated 
environments  there  need  to  be  simulated  with  a method  of  analysis  that  1)  uses  a valid  intra-room 
fire  dynamics  modeling  approach  and  2)  can  be  implemented  in  the  compartment  fire  model  to  study 
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fire  scenarios  in  facilities  which  include  other  "standard,"  i.e.,  small-to-moderate-I/^/,  room  elements. 
The  purpose  of  this  paper  is  to  develop  a framework  for  such  an  analysis. 

Once  a method  of  analyzing  large- I/c/  spaces  is  developed  and  verified,  it  could  be  incorporated  into 
an  existing  multi-room  fire  model.  The  revised  model  would  then  be  capable  of  describing  smoke 
movement  through  shaft-like  spaces  that  may  be  included  in  the  design  of  a particular  facility  of 
interest.  The  strategy  for  carrying  out  simulations  with  the  revised  fire  model  would  involve 
categorizing  each  room  of  a modeled  facility  as  either  a small-to-moderate-I/^/  space  or  a large- I/i/ 
space,  where  results  of  future  research  are  expected  to  lead  to  appropriate  rules  for  such 
categorization.  Then,  when  carrying  out  a simulation,  the  fire  model  would  invoke  the  traditional, 
two-layer,  isolated-plume,  etc.  model  equations  for  the  small-to-moderate-I/^f  spaces  and  the  new 
model  equations  for  large- I/J  spaces. 

Note  that  an  even  more  advanced  strategy  might  consider  an  additional,  third  class  of  facility  space, 
designated  as  a "moderate- I/c?"  space.  The  analysis  of  the  fire  environment  in  the  "moderate- 
space  would  be  designed  to  bridge  the  gap  between  small-  and  large- I/c?  spaces  by  including 
important  modeling  features  of  both. 

COMBINED  BUOYANCY-  AND  VENTILATION-DRIVEN  FLOW  OF  A PERFECT  GAS 
THROUGH  A LONG  VERTICAL  SHAFT 

Consider  flow  of  a perfect  gas  through  a vertical  shaft  of  length  L.  Let  the  section  of  the  shaft  have 
a characteristic  dimension,  d,  and  assume  that  the  shaft  is  long  in  the  sense  that  Ljd  > > 1 . p and 
r,  the  density  and  temperature  of  the  gas  averaged  across  the  section,  vary  along  the  length  of  the 
shaft.  In  general,  the  shaft  is  ventilated;  inflow  or  outflow  of  the  gas  can  occur  at  vents,  in  the  walls 
of  the  shaft,  or  at  its  ends.  Inflowing  gas  has  specified  density  and  temperature,  and  Ty^, 
respectively,  which  can  vary  along  the  shaft.  Outflowing  gas  is  assumed  to  be  at  density  and  tempera- 
ture, p and  r,  respectively.  The  local  rate  of  mass  addition  due  to  ventilation  is  (rate  of  mass 
added  to  the  shaft  per  unit  volume).  Heat  is  transferred  between  the  gas  and  the  shaft  surface  at  the 

rate  (rate  of  heat  transferred  to  the  gas  in  the  shaft  per  unit  volume). 

Assume  that  the  only  possible  significant  component  of  flow  velocity  is  along  the  axis  of  the  shaft  in 
the  vertical  z direction  and  let  its  average  across  a section  be  V^.  Then,  conservation  of  mass  and 
energy,  and  the  equation  of  state  leads  to: 


3 


0p/5f  + 3(K^p)/3z  - 


(1) 


pC^(dT/dt  ^V^dTldz)  = - d(kdT/dz)ldz  ^ CJiTy^  - T)m'^  + q'^  (2) 


= Pamb’^amb  = constant 


(3) 


where  is  the  specific  heat  at  constant  pressure,  k is  the  thermal  conductivity,  and  p^^  and 
are  the  density  and  temperature  at  some  ambient  reference  state. 

Note  that  the  formulation  of  Eq.  (3)  uses  the  uniformly  valid  assumption  that  variations  in  pressure 
throughout  the  length  of  the  shaft  are  negligible  compared  to  the  characteristic  pressure  of  the 
ambient  reference  state.  At  a point  in  the  analysis  where  momentum  considerations  are  required, 
for  example  in  the  computation  of  pressure  differences  across  vents  for  the  purpose  of  estimating 
vent  flow  rates,  estimates  of  the  relatively-small  shaft  pressure  variations  (e.g.,  the  contribution  of 
hydrostatic  pressure  variation)  would  be  important  and  they  must  be  estimated  properly. 

Using  Eqs.  (1)  and  (3)  in  Eq.  (2)  leads  to  the  following  simplified  alternative  to  Eq.  (2): 


dVJdz  = [-  d(kdTldz)ldz  - C;ryg^m'^  - ambT' amb^ p) 


(4) 


TURBULENT  FLUCTUATIONS  IN  THE  SHAFT  FLOW 
The  'Hirbulence  Equations 

Eqs.  (1),  (3),  and  (4)  are  now  revised  to  account  for  the  effect  of  possible  turbulent  fluctuations  in 
the  shaft  flow.  Time-averaged  ("barred")  and  fluctuating  ("primed")  components  of  the  variables  are 
introduced 

p = p + p';  T =T  * T\  * V'  (5) 


Reynolds  averaging,  neglect  of  the  generally  insignificant  effect  of  heat  conduction  through  the  gas, 
along  the  length  of  the  shaft,  as  derived  from  the  time-averaged  temperature  distribution,  and  the 
assumption  that  time-averaged  values  of  p and  T satisfy  the  equation  of  state  lead  to  the  following 
modified  equation  set: 
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ap/ar  + v^dpjdi  * aCVjp')/*  = m^(l  - p/Pke^)  - (PlPAm'^imKI'AMB^p) 


(6) 


dVJdz  - (S^^YEtij^vEtfr  ■*■  ^ht^I^Pamb^amb^p^ 


= Pamb'^aub  = constant 


(7) 

(8) 


Turbulent  fluctuations  leading  to  significant  values  of  the  term  in  Eq.  (6)  will  occur  along  the 
length  of  a vertical  shaft  where  there  are  buoyancy-generated  instabilities  in  the  vertical  density 

distribution,  i.e.,  where,  locally,  the  density  is  increasing  with  elevation,  dpfdz  > 0.  Indeed,  such 
instabilities  are  the  driving  force  for  the  turbulent-like  mixing  phenomenon  which  is  of  particular 

interest  here.  Also,  shaft  flow  scenarios  of  present  interest  are  such  that  is  small  enough  as  never 

to  lead  to  turbulent  fluctuation  enhancements  that  would  significantly  affect  the  value  of  For 
example,  at  elevations  along  the  shaft  where  the  gas  is  stably  stratified,  i.e.,  where  dp/dz  ^ 0,  there 

will  be  no  buoyancy-driven  turbulence  and  it  is  reasonable  to  assume  that  3(F/pO/dz  in  Eq.  (6)  can 
be  neglected  completely. 

Solutions  to  Eqs.  (6)-(8)  with  appropriate  initial  and  boundary  conditions  would  generally  provide  a 
description  of  the  fire  environment  that  develops  in  ventilated,  shaft-like,  room  configurations. 
However,  to  actually  implement  this  equation  set  it  is  necessary  to  develop  supplementary  model 

equations  to  represent  F/p'  of  Eq.  (6)  and  the  source  terms,  and  of  Eqs.  (6)  and  (7). 

The  focus  of  the  remainder  of  the  present  work  is  on  obtaining  a general  representation  for  F^p^ 

Experiments  in  Shafts  With  Unstable  Density  Configurations;  A Qualitative  Result 

During  experiments  in  vertical  tubes  of  diameter  d,  involving  initially  unstable  configurations  of  salt- 
water over  fresh-water  and  prior  to  any  end  effects.  Cannon  and  Zukoski  [1]  have  observed  fluid 
mixing  by  buoyancy-generated  instabilities.  This  involved  overturning  and  fluctuating  eddies  of 
characteristic  dimension  d.  Also,  as  observed  in  similar  unstable  systems  by  Cannon  and  Zukoski  [1] 
(heavy  gas  over  light  gas),  Baird  and  Rice  [2]  (rising  bubble  columns),  and  Gardner  [3]  (salt- 
water/fresh-water  systems  and  experiments  on  water  rising  through  more  dense  carbon  tetrachloride), 
the  propagation  of  buoyancy-driven  turbulent  fluctuations  through  the  buoyantly  unstable  column 

behaves  phenomenologically  as  a diffusion-like  process.  This  is  represented  by  the  d(V^p^/dz  term 
in  Eq.  (6). 
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A Representation  for 

Shaft  flow  systems  involving  negligible  V^.  Compressibility  was  negligible  in  all  the  above-referenced 
experiments.  Therefore,  was  negligible  everywhere,  and  it  would  appear  that  insights  on  a 

practical  and  general  representation  for  in  Eq.  (6)  can  be  obtained  from  the  results  of  these 
relatively-simple  negligible-  systems. 

Consider  shaft  flow  systems  of  Eqs.  (6)-(8)  with  negligible  and  negligible  mass  addition  or  heat 
transfer  along  the  length  of  the  shaft.  For  such  systems,  Eq.  (6)  becomes 

ap/ar  + a(K,p^/az  = o (^) 


Defining  D as  the  eddy  diffusivity  for  buoyancy-driven  turbulence,  Eq.  (9),  with 


f'p'  = - oap/az 


(10) 


has  been  used  in  each  of  the  above-referenced  systems.  It  is  also  used  by  Cannon  and  Zukoski  [1] 
in  the  modeling  of  a perfect  gas  air  system,  i.e.,  low  temperature  over  high-temperature  air,  with  non- 
zero q'!^.  In  Eq.  (10),  D is  expected  to  be  an  increasing  function  of  0p/5z  for  dp/5z  > 0,  and  to  be 

zero  for  8p/3z  ^ 0 (i.e.,  buoyancy-driven  turbulent  diffusion  is  suppressed  at  elevations  where  the 
fluid  column  is  stably  stratified). 

Cannon  and  Zukoski  [1]  expect  Z)  to  be  proportional  to  a characteristic  velocity  of  the  turbulent 
fluctuations,  and  to  satisfy 


D - w'l  = [(g/p)ap/az]''^/^ 

where  / is  a characteristic  length  of  the  fluid  mixing  process.  Here,  / is  taken  to  be  the  characte- 
ristic dimension  of  the  observed  buoyancy-driven  eddies,  namely,  d. 

^ _ I Kd\(glp)dpldz\'^  if  dpjdz  > 0 
1 0 i/  dpjdz  ^ 0 
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where  is  a constant  and  where  the  form  of  Eq.  (11)  is  generally  equivalent  to  that  presented  by 
Baird  and  Rice  [2]  and  Gardner  [3]. 


A determination  of  ^ ; a steady-state  experiment.  To  determine  K,  Gardner  [3]  considered  the 
particular  problem  of  quasi-steady  turbulent  diffusion  in  a long  shaft  connecting  two  relatively  large 
vessels  of  incompressible  fluids,  with  density  of  the  fluid  in  the  top  and  bottom  vessels,  and 
Pbot"^  Pro,,  respectively;  where  a p = pro,  - Pjorissuchthat  Ap/(pro,  + Psot^  <<  1-  To  insure 
that  = 0,  at  least  one  of  the  vessels  is  fully  enclosed.  This  problem  was  studied  experimentally 
by  Mercer  and  Thompson  [4]  with  salt-water/freshwater  systems  and  circular  shafts  of  diameter,  d. 

For  the  above  class  of  problem,  p « = (pj^^  + 

D = Kd^[(sAplp^y^)dc^^dz]''^;  - DdCgJdz  = Ddc^^dz  = D/L  = 4(?/(rt<i^)  (12) 


where  Q is  the  volume  flow  rate  of  the  fluid  in  the  lower  vessel  diffusing  up  through  the  shaft  and 
into  the  top  vessel  and  where  Cj^p,  are  the  concentrations  in  the  shaft  of  the  top  and  bottom 
fluids,  respectively. 

With  quasi-steady  systems  of  salt  water  over  fresh  water,  Mercer  and  Thompson  [4]  determined^ 
of  fresh  water  from  the  lower  vessel  by  measuring  salt  concentrations  and  the  rates  of  mass  decrease 
of  the  salt  water  in  a top  vessel.  Using  these  results  in  Eqs.  (12),  Gardner  [3]  determined 

K ^ = 0.56  (13) 


In  contrast  to  this,  Baird  and  Rice  [2]  analyzed  the  literature  on  rising  bubble  columns  and  related 
systems  where  gas  flow  rates  are  large  enough  for  the  bubble  wakes  to  interact  strongly  and  to 
generate  turbulence  in  the  opposing  fluid  motion.  For  such  systems  they  determined  the  "order-of- 
magnitude  estimate” 


K = Kpp  = 0.21 


(14) 


UNSTEADY  EXPERIMENTS  TO  VERIFY  EQS.  (9)-(ll)  AND  THE  VALUE  OF  K 
Description  of  the  T\vo  Phase  Experiments 

The  unsteady  experiments  of  Cannon  and  Zukoski  [1]  included  the  configuration  of  Figure  2 
involving  a tube  of  length  I,  closed  at  the  upper  end,  z = L,  and  inserted  into  a relatively  large 


7 


vessel  of  density  The  tube  was  initially  filled  with  fluid  of  density  The 

experiments  were  initiated  by  allowing  the  fluid  in  the  tube  to  mk  with  the  fluid  in  the  lower  vessel. 
The  mixing  phenomena  in  the  tube  involves  two  phases,  Phases  I and  II.  Phase  I begins  when  the 
experiment  is  initiated.  A diffusion  front  at  elevation  z = ZppQi^(t),  = 0)  =0,  rises  in  the 

tube,  and  in  0 ^ z ^ there  is  mixing  due  to  non-zero  dpjdz.  Phase  I is  completed  and 

Phase  II  is  initiated  at  time  r = when  Zp^Qj^j-  = L,  i.e.,  at  the  time  of  arrival  at  the  top  of  the 

tube  of  the  diffusion  front,  when  the  effects  of  the  diffusive  mixing  begin  to  modify  significantly  the 
density  at  the  upper  end  of  the  tube.  During  the  remainder  of  Phase  II,  for  t > tpj^j^,  the  unstable 
density  distribution  drives  turbulent-like  mbdng  eddies  through  the  length  of  the  tube.  The 
distribution  eventually  relaxes  to  the  uniform  state,  p = 

Data  from  both  Phases  I and  II  of  the  Cannon  and  Zukoski  experiments  [1]  can  be  used  to  verify  the 
model  Eqs.  (9)-(ll)  and  the  value  of  K.  The  next  section  addresses  the  Phase  I phenomena. 

The  Initial  Value  Problem  for  the  Diffusion  Front:  Phase  I 

Problem  Formulation  and  Similarity.  From  Eqs.  (9)-(ll),  the  initial  value  problem  for  Phase  I of 
the  above-referenced  experiments  is 


d6  _ d [ 1 I SB  \3/2j 


(15) 


6(C  0,  T > 0;  €)  0;  1)  ^ front’  ^ ^ front  ® 


(16) 


e ^ ~ ^ Klegld]^t;  € ^ 4p/p,or 

(PjOF  - Pbot) 


(17) 


Note  that  Eqs.  (15)-(17)  specify  p(z  = 0)  = PpQp.  However,  as  discussed  by  Cannon  and  Zukoski 
[1],  p at  the  base  of  the  tube  will  always  be  somewhat  greater  than  p^^,  i.e., 

0(C  = 0;  T > 0)  = 6(t)  where  6 presumably  satisfies  0 < 6 <<  1.  Thus,  the  indicated  boundary 
condition,  0(C  = 0;  t > 0)  = 0,  of  Eq.  (16)  is  only  approximately  correct.  In  particular  the  solution 
to  Eqs.  (15)-(17)  can  be  expected  to  provide  a good  simulation  of  experimental  values  only  for 
moderate-to-large  values  of  C and  CppoNT' 

Introducing  the  similarity  variable  q , it  can  be  shown  that 
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0 = 0(r|;  €);  where  t|  = 


(18) 


and  that  the  problem  of  Eqs.  (15)-(17)  can  be  reduced  to 


T1 


dr\ 


5 d 

1 

2 dx] 

[(1  + eO)*^ 

\dr]  / . 

; 0 ^ 0 ^ 1; 


6(11  = 0)  = 0;  -^(6  = 1) 
dy\ 


'^(’1  “ '^FROSt)  ~ '^PRONT  ~ '^FROtrM^ 


(19) 


where  is  the  smallest  value  of  r\  satisfying  the  above  condition. 

Once  ^ppoffpi^)  is  obtained  in  the  solution  of  Eqs.  (19),  the  time-dependent  position  of  the  diffusion 
front  can  be  determined  from 

<^FROm(f’  «)  = ^FRONl(^>  6)/</  = 


Solution  for  small  e : verification  with  salt-water/fresh-water  experiments.  For  the  salt-water/fresh- 
water experiments  of  Cannon  and  Zukoski  [1],  e always  satisfies  € < 0.2.  Therefore,  thee  = 0 
solution  to  Eqs.  (15)-(17)  should  provide  a good  simulation.  This  solution  is  found  to  be 


0(TT,  € = 0)  = if— H— If— il— 1 - ilf_Jl_j 

8 V Tj  ppQJsfJ'  ^ ^o/’wrr'  3 V T1  co/^Rrr' 


front'  ^ FRONT 


+ 5 


, 0 ^ 


FRONT 


^ 1; 


(21) 


\ FRONT 


= ’lF*D«T(e  = 0)  = 15''^/8‘'5  = 3.3499... 


0 < € « 1)  = Zj,«o^(f;  0 < e « D/d  = (15^'ys'>^)[K(€gld)'>^tf^  (22) 


Eq.  (22)  is  plotted  in  Figure  3 for  the  K's  of  Eqs.  (13)  and  (14).  (The  K = 0.438  plot  will  be 
explained  below  in  the  discussion  of  Phase  II.) 

The  salt-water/fresh  water  experiments  of  Cannon  and  Zukoski  [1]  cover  the  d,  Ljd  range 
0.030m  ^ d ^ 0.181m,  7.6  ^ Lfd  ^ 40.3.  For  many,  but  not  all  of  these  experiments  they  tabulate 
(Table  6-2  of  [1])  and/or  plot  (Figure  6.13  of  [1])  These  data  are  included  in  Figure  3.  In 
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the  figure,  note  that  when  K is  taken  to  be  Kq  = 0.56  the  larger  values  of  the  Cprqnt  well- 

predicted  by  Eq.  (22).  (Recall  that  this  value  of  K was  also  deduced  from  experiments  in  a salt- 
water/fresh-water  system.)  It  is  also  noteworthy  that  data  for  the  largest  Ljd  experiments 

Ljd  = 29.3  and  40.3,  are  not  reported  in  [1]  and,  therefore,  are  not  included  in  Figure  3.  However, 
data  from  these  experiments  up  to  Cpront  ~ reported.  For  times  prior  to  the 

K = 0.56  plot  of  Figure  3 continues  to  predict  well  this  data.  Also,  because  of  the  aforementioned 
approximate  nature  of  the  q =0  boundary  condition  of  Fqs.  (19),  and  because  the  characteristic  d- 
dimension  of  the  turbulent-like  eddies  are  associated  with  C intervals  of  the  order  of  1,  small-to- 
moderate  values  of  expected  to  be  well-simulated  by  the  Fq.  (22)  solution. 

Solution  for  moderate-to-larae  e : verification  with  heaw-gas/]i2ht-2as  experiments.  Simulations  of 
the  unsteady  heavy-gas/light-gas  experiments  of  Cannon  and  Zukoski  [1]  require  solutions  to  Fqs. 

(19)  with  moderate-to-large  values  of  e.  To  study  such  cases  it  is  convenient  to  solve  the  following 
alternative  coupled  set  of  equations  for  0 , \|; , and  4) , 


d^ 


4^ 

15/ 


+ €0)^ 


13^(1  + €0) 


0^0 

dx]  dr\ 


^ 1 


where  0(r|  = 0)  = i|;(r)  = 0)  = 0 


and  where  r(€)  = (|)(t| 


0 = 


= 0)  = — (q  = 0)  is  determined  such  that 
dr\ 

= 0 at  some  minimum  q = 


(23) 


Solutions  for 


4^(ti;  e)  = p0(G);  e)d(ji 

J 0 


(24) 


will  be  used  below  in  the  Phase  II  analysis. 

Numerical  integration  of  Fqs.  (23)  was  carried  out.  This  yielded  0(q;  e),  r(€),  and 

ij;(q;  e)  for  0 < € ^40.  Some  results  are  plotted  in  Figures  4a,  4b,  and  4c,  respectively.  Figure  4a 
also  includes  a plot  of  the  0(q;  e = 0)  result  of  Fq.  (21).  As  can  be  seen,  although  there  is  a 
significant  variation  in  ^ppour^^^’  are  relatively  insensitive  to  differences  in 

€.  plotted  in  Figure  4d. 
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Cannon  and  Zukoski  [1]  report  on  Figure  2-type  heavy-gas/light-gas  experiments  with  d = 0.152»i, 
LJd  = 12  and  with  € = 0.22,  0.38,  and  4.0.  Figure  6.14  of  [1]  presents  Phase  II  data  on 
0(C  = Lfdf  T > ^xjijijvAL  » However,  Phase  I data  for  these  experiments  are  not  provided. 

Also,  note  that  estimates  of  by  extrapolation  of  the  Phase  II  data  to  0 = 1 are  subject  to 

significant  error  because  increases  in  e lead  to  decreases  in  ^arrival  ^ climated  from 

Eq.  (26)  to  be  of  the  order  of  5s  when  e = 4.0],  and  because  of  the  "slow  response  time  of  the 
system  used  to  measure  the  density  of  the  gas  mixtures  [1]." 

Relaxation  of  the  Diffusion  Front:  Phase  II 

The  initial  value  problem  for  Phase  II.  Subsequent  to  r = during  the  relaxation  of  the 

diffusion  front,  0 is  no  longer  characterized  by  the  similarity  variable.  Rather,  it  must  be  determined 
from  a solution  of  Eq.  (15)  subject  to 


e(C  = 0,  T > e)  = |?(f  = w T > e)  = 0 


ac 


e(C,  t €)  - z(C;  €)  = e(n  = €) 


(25) 


^ARRIVAL 


(26) 


An  approximate  solution.  Partial  verification  of  Eqs.  (9)-(ll)  during  Phase  II  can  be  obtained  by 
an  integral-type  approximate  solution  to  Eqs.  (15)  and  (25).  Toward  this  end,  Eq.  (15)  is  first 
integrated  in  C in  the  interval  0 ^ C ^ Lfd.  Using  the  boundary  conditions  of  Eqs.  (25),  this  leads 
to 


d_ 

dx 


r^0(C,  t;  e)dC 


Vacic 


,3/2 


(27) 


Now  choose  the  following  approximate  solution  form,  which  satisfies  exactly  the  initial/boundary 
conditions  of  Eq.  (25). 

0 - 0(C>  = ^ARRIVAL*  “ ^ARRIVAL*  ~ ^ (^) 


Using  Eq.  (28)  in  Eq.  (27)  eventually  leads  to 
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TTt;  £.  Ud)  = ([0(e)/(LW^'^(T 

‘^ARRIVAl)  U ^ 


(29) 


where  Q(€)  is  plotted  in  Figure  5. 

Solution  for  small  e ; salt-water/fresh-water  experiments.  For  small  e,  Eqs.  (21),  (22),  and  (29)  lead 
to 


Q(€  = 0)  = {l5^^/2^'^)/n  = 1.8672... 


(30) 


Using  Eqs.  (28)-(30),  the  Phase  II  solution  at  z = L for  e ^ 0 and  Ljd  = 15.3  is  plotted  in  Figure 
6.  TTiis  is  the  approximate  solution  that  would  be  used  to  simulate  Figure  2-type  salt-water/fresh- 
water  experiments  of  Cannon  and  Zukoski  [1].  The  most  accessible  Phase  II  data  from  these 
experiments  are  presented  in  Figure  6.9  of  [1]  for  the  conditions:  d = 0.089m,  Ljd  = 15.3, 
e = 0.050.  Reducing  these  data  for  plotting  in  Figure  6 requires  the  value  of  K.  Two  values  were 
used,  and  corresponding  reduced  data  sets  were  obtained  and  plotted.  The  first  is  X"  = = 0.56. 

The  second  is  that  value  of  K which  provides  the  least  squares  fit  between  the  data  and  the 
theoretical  solution  curve.  The  latter  value  was  found  to  be  A = 0.438 , and,  as  can  be  seen  in  the 
figure,  this  provides  an  excellent  match  between  the  theory  and  the  experiment  during  Phase  II. 

Using  the  "best-fit"  K = 0.438,  a plot  of  the  corresponding,  theoretical.  Phase  I-solution  has  also 
been  included  in  Figure  3.  As  can  be  seen,  this  value  of  K and  the  corresponding  solution  provides 
a somewhat  better  match  to  the  Phase  I data  than  does  that  of  the  K = 0.56  solution. 

Solutions  for  moderate  e;  heaw-gas/light-oas  experiments.  Using  Eqs.  (28)-(29),  the  Phase  II 
solutions  at  z = L for  € = 0.22,  0.38,  and  4.0  are  plotted  in  Figure  7.  These  approximate  solutions 
are  to  be  compared  to  the  heavy-gas/light-gas  data  of  Figure  6.14  of  Cannon  and  Zukoski  [1].  Again, 
reduction  of  the  data  for  presentation  in  this  figure  requires  a value  of  K.  This  is  chosen  to  the 
previously  determined  "best-fit"  value,  K = 0.438 . As  can  be  seen,  the  comparison  between  theory 
and  experiment  is  good,  especially  for  € = 0.22  and  0.38. 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

A limitation  of  existing  zone-type  compartment  fire  models  is  the  lack  of  a satisfactory  means  of 
simulating  time-dependent  fire-generated  environments  in  ventilated  shaft-like  spaces  with  large 
height-to-span  {Ljd)  ratios.  For  such  configurations,  the  traditional  time-dependent,  stably-stratified, 
spatially-uniform,  two-layer  description  of  a room  fire  environment,  involving  relatively-small-volume 
quasi-steady  buoyant  plumes  is  not  valid.  Rather,  experimental  observations  previously  reported  in 
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the  literature  suggest  a time-  and  elevation-dependent  density  distribution,  where  any  intervals  of 
instability  (i.e.,  density  increasing  with  elevation)  tend  to  be  removed  by  the  action  of  overturning  and 
fluctuating  turbulent-like  mixing  eddies  of  characteristic  dimension,  d.  It  has  also  been  observed  that 
on  the  macroscopic  level,  the  mbdng  phenomenon  has  the  features  of  a diffusion-like  process. 

Consistent  with  the  above  ideas  and  observations,  a set  of  model  equations  to  simulate  combined 
buoyancy-  and  ventilation-driven  flow  of  a perfect  gas  through  a long  vertical  shaft  was  developed  and 
presented  as  Eqs.  (6)-(8),  (lO)-(ll).  Required  to  implement  and  solve  this  equation  set  is  the  value 
of  a constant  K.  This  is  the  coefficient  of  the  equations’  turbulent  diffusivity  term  that  leads  to 
modifications  in  the  density  distribution  in  any  regions  along  the  length  of  the  shaft  where  an  unstable 
density  distribution  prevails. 

The  model  equations  were  studied  for  a special  class  of  problem  involving  no  net  ventilation  flow 
velocity  and  no  wall-to-gas  heat  transfer.  Under  such  conditions  all  changes  in  density  distribution 
are  a direct  result  of  the  above  diffusion-like  mixing  mechanism.  Also,  the  equations  are  directly 
applicable  to  systems  involving  two-component  fluids  or  gases. 

For  the  above  class  of  problem  it  was  determined  that  with  K = 0.44 , data  from  relevant,  unsteady. 
Figure  2-type  salt-water/fresh-water  experiments  of  Cannon  and  Zukoski  [1]  are  well-simulated  by 
the  equation  set,  Eqs.  (9)-(ll).  For  Ljd  ^ 15,  excellent  correlation  between  theory  and  experiment 
was  found  for  both  Phases  I and  II  of  the  experiment.  (Phase  I is  when  the  diffusion  front,  initiated 
at  the  bottom,  moves  up  the  shaft  through  the  initially  undisturbed  liquid.  Phase  II  begins  when  the 
front  reaches  the  top  of  the  closed  shaft  and  continues  until  the  unstable,  top-to-bottom,  density 
distribution  relaxes  asymptotically  in  time  to  a uniform  state.)  Data  from  additional  Figure  2-type 
experiments  involving  heavy-gas/light-gas  systems  and  relatively-large,  initial,  density  differences 
(€  = 0.22,  0.38,  and  4.0)  were  also  well-simulated  by  the  Eq.  (9)-(ll),  K = 0.44  equation  set. 

The  above  results  lead  to  some  confidence  in  the  validity  of  describing  the  diffusion-like  mixing 
processes  of  interest  by  Eqs.  (10)  and  (11)  with  K = 0.44.  In  order  to  form  an  even  greater  level 
of  confidence,  additional  comparisons  between  theory  and  experiment  are  necessary,  especially  in 
configurations  different  from  that  of  Figure  2. 

In  view  of  the  above  it  is  recommended  that  the  full  equation  set,  Eqs.  (6)-(8),  (lO)-(ll),  with 
K = 0.44  be  used  in  room  fire  models  to  describe  the  development  of  fire  environments  in  ventilated, 
shaft-like,  room  configurations.  This  will  require  further  development  of  the  heat  transfer  elements 

of  Eqs.  (6)  and  (7),  i.e.,  a verified  method  of  simulating  the  term.  The  new  large- I/t/  model 
equations  will  be  subject  to  verification  with  data  from  experiments  on  hot-air/cold-air  systems  and 
on  fire-driven  smoke  flows,  both  of  which  involve  gas-to-surface  heat  transfer  exchanges.  Such 
verification  is  required  before  the  new  shaft  model  can  be  used  with  confidence  in  professional 
practice. 
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NOMENCLATURE 


specific  heat  at  constant  pressure 

^BOT^  ^TOP 

concentrations  of  top,  bottom  fluids 

D 

eddy  diffusivity 

d 

characteristic  span  of  a room  or  section  of  a shaft,  e.g.,  diameter 

8 

acceleration  of  gravity 

K 

constant,  Eq.  (11) 

^G’  ^BR 

Eqs.  (13),  (14) 

k 

thermal  conductivity 

L 

characteristic  height  of  a compartment  space,  e.g.,  length  of  a shaft 

1 

characteristic  length  of  the  fluid  mixing  process 

. /// 

^VENT 

rate  of  mass  added  to  the  shaft  per  unit  volume 

Q 

volume  flow  rate 

. /// 

^HT 

rate  of  heat  transferred  to  the  gas  in  the  shaft  per  unit  volume 

T 

temperature  in  the  shaft  averaged  across  the  section 

T 

^VENT 

T of  gas  flowing  through  a vent  and  into  the  shaft 

j,  j.f 

time-averaged,  fluctuating  component  of  T 

t 

time 

^ARRIVAL 

t when  ZpjfQffj  = L 

K 

flow  velocity  in  z direction  average  across  a section 

K K' 

time-averaged,  fluctuating  component  of 
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w' 

characteristic  velocity  of  the  turbulent  fluctuations 

z 

coordinate  in  upward  direction 

^FROST 

z of  diffusion  front 

r 

Eq.  (23) 

Ap 

Prop  “ Ppor 

£ 

Eq.  (17) 

c 

Eq.  (17) 

^ FRONT 

C at  the  diffusion  front 

Eq.  (18) 

^ FRONT 

q at  C = CpRONT 

e 

Eq.  (17) 

P 

density  in  the  shaft  averaged  across  the  section 

pAi/B’  Pfejvt 

p of  ambient,  gas  flowing  through  a vent  and  into  the  shaft 

P»  P^ 

time-averaged,  fluctuating  component  of  p 

Pafe 

(Prop  PpoT^/^ 

Ptop»  Psor 

p of  fluid  in  top,  bottom  of  the  vessels;  initial  p at  top,  bottom  of  shaft 

T 

Eq.  (17) 

“^ARRIVAL 

Eq.  (26) 

i> 

Eq.  (23) 

Eqs.  (23),  (24) 

Q 

Eq.  (29) 

(0 

dummy  variable,  Eq.  (24) 
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Figure  1.  Sketch  of  a facility  where  traditional,  two-layer,  zone-type,  modeling  concepts 
are  not  uniformly  applicable. 
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t = 0 


0<t<t 


ARRIVAL 


Figure  2.  Sketch  of  the  simulated  experiments  of  Cannon  and  Zukoski  [1], 
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10  20  100  200 


[(g/d)(Ap/pBOT)]'''t 


Figure  3.  Position  of  the  diffusion  front  in  Phase  I:  small-  e solution  of  Exj.  (22)  for  different 
values  of  and  values  measured  in  the  salt-water/fresh-water  experiments  of 
Cannon  and  Zukoski  [1]. 
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FRONT 


(a) 


Figure  4.  Phase  I solutions  of  Eqs.  (23),  including  moderate-to-large  e:  (a)  0(ti;  c); 
(b)  r(e);  (c)  and  (d)  e]. 
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(b) 


Figure  4.  Phase  I solutions  of  Eqs.  (23),  including  moderate-to-large  e:  (a)  0(t|;  e); 
(Cont’d)  (b)  r(€);  (c)  and  (d)  thraowC®); 
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Figure  4.  Phase  I solutions  of  Eqs.  (23),  including  moderate-to-large  e:  (a)  0(ti;  e); 
(Cont’d)  (b)  r(e);  (c)  and  (d)  e],  " 
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Figure  4.  Phase  I solutions  of  Eqs.  (23),  including  moderate-to-large  e:  (a)  0(t|;  e); 
(Confd)  (b)  r(€);  (c)  and  (d)  ^1- 
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Figure  5.  Plot  of  Q(€)  of  Eq.  (29)  from  Phase  II  solutions  of  Eqs.  (27)  and  (28). 
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0(^=L/d,  t:e  = 0,  L/d  = 15.3) 


Figure  6. 


■ "^ARRIVAL 


Plots  for  Phase  II:  Eqs.  (28)-(29)  at  z = L for  € =0;  salt-water/fresh-water 
data  of  Cannon  and  Zukoski  [1]  using  K = 0.56  and  least-squares-fit  value, 
K = 0.438. 
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0(^  = L/d,  t;  8,  L/d  = 12) 
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Figure  7. 
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Plots  for  Phase  II:  Eqs.  (28)-(29)  at  z = L for  € = 0.22,  0.38,  and  4.0; 
heavy-gasAight-gas  data  of  Cannon  and  Zukoski  [1]  for  e = 0.22,  0.38,  and 4.0 
using  K = 0.438. 
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